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Abstract 

The spectrum of the massive Schwinger model in the strong coupling region 
is obtained by using the light-front Tamm-Dancoff (LFTD) approximation 
up to including six-body states. We numerically confirm that the two-meson 
bound state has a negligibly small six-body component. Emphasis is on the 
usefulness of the information about states (wave functions) . It is used for iden- 
tifying the three-meson bound state among the states below the three-meson 
threshold. We also show that the two-meson bound state is well described by 

the wave function of the relative motion. 
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I. INTRODUCTION 



In a previous paper [JJ, we investigated the massive Schwinger model [|2],|3| with SU(2)f 
in the light-front Tamm-Dancoff (LFTD) approximation [|]|| up to including four-body 
states. We showed, by examining the wave functions, that the lightest isosinglet state can 
be regarded as a bound state of two "pions." This observation naturally led us to the answer 
to the question raised by Coleman || why it is so light. The LFTD approximation has been 
proved to be one of the most powerful non-perturbative methods to date in the investigation 
of relativistic bound states, at least in two dimensions, although we have to face the difficult 
renormalization problem in higher dimensions. 

The validity of the LFTD approximation is based on the plausible hope that the 
sea quark/gluon contributions are small in the light-cone quantization because pair cre- 
ations/annihilations are suppressed |J. Typically, the lightest particles are expected to be 
in the valence states. It is generally true in the models so far investigated. The above 
mentioned state (the bound state of two "pions") is an important exception. With this 
exception, one might think that such a state would have non-negligible many-body compo- 
nents too. It is one of our purposes of this paper to show numerically that it is unlikely by 
examining the single-flavor model. 

We also investigate the three-meson bound state of the single-flavor model. Its existence 
has been discussed by Coleman |§ by using the bosonization technique. He showed that, 
in the strong coupling limit with the zero vacuum angle, there exists a stable three-meson 
bound state and it is unstable when the vacuum angle is non-zero. We look for a candidate 
which can be interpreted as a bound state of three mesons in our numerical results. 

In order to investigate these problems, it is necessary to do LFTD calculations up to 
including six-body states. Such calculations are very hard without any technical refinement. 
To make these calculations feasible we have made two points: (1) We take a simple set of 
basis functions in order to reduce CPU time. A clever choice of basis functions is essential 
as in quantum chemistry calculations. Note that our choice of basis functions in a finite 
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domain will be also useful for higher dimensions. Even in higher dimensions the longitudinal 
momenta pf of constituents are restricted to a finite domain < pf < P + with Y^iPt — P + , 
where P + is the total momentum. (2) The three-meson bound state, if it exists, must be in 
the continuum unless it is lighter than two mesons, and is therefore apparently difficult to 
find. We can however find a candidate among several states by looking at the wave functions. 
The points are that a three-meson state must be charge conjugation odd and that below the 
three-meson threshold, six-body components should be very small except for three-meson 



bound states. A more detailed discussion is given in Sec. |Tj. 

We emphasize that the information about states (wave functions) is very useful. It is 
used for identifying the three-meson bound state, as is said above. As another example, 
we introduce the wave function of the relative motion of the two-meson bound state and 
try to describe the bound state in terms of the wave function. Although the concept of 
"relative motion" of a relativistic bound state is somewhat awkward, we however find that 
the two-meson bound state is well described in terms of the wave function of the relative 
motion, in the sense that a smaller set of basis functions motivated by the concept of the 
relative motion gives a good approximation. It gives us a qualitative picture of the bound 
state. 

We summarize the results: (1) The masses of the lowest states do not change even if we 
include six-body states. (2) In particular, the state which can be regarded as a bound state 
of two mesons has a negligible six-body component. (3) We find a candidate for the bound 
state of three mesons. (4) The wave function of the relative motion of the two-meson bound 
state describes the bound state well. We can have a picture that in the strong coupling 
region it is loosely bound, while in the weak coupling region it is tightly bound, compare to 
the size of the meson. 

The massive Schwinger model J7|,p| has been discussed by many authors in the light-cone 
quantization. Bergknoff || did the first LFTD calculations. Mo and Perry [jXOj] refined his 
calculations by the use of basis functions. Their calculations include only up to four-body 
states. Eller, Pauli, and Brodsky [Tl]|l2j discussed the massless and the massive Schwinger 



models in the discretized light cone quantization (DLCQ). Our work is based especially on 
the papers by Bergknoff, and Mo and Perry. We try to keep our notation as close as possible 
to that of our previous paper [I]. 

In Sec. [TI|, we present basic facts and formulas on the massive Schwinger model to make 
this paper self-contained. The model is quantized on the light cone and the Tamm-Dancoff 
truncation is made up to including six-body states. The wave functions are expanded in 
terms of a new set of basis functions. The numerical results are shown in Sec. pT[ We identify 
two-meson and three-meson bound states. The two-meson bound state is shown to have a 
negligibly small six-body component. The three- meson bound state is charge conjugation 
odd and has a large six-body component compare to those of other states below the three- 
meson threshold. In Sec. [TV], we introduce a meson operator which (approximately) creates 
a meson from the vacuum. By using the meson operator, we also introduce the wave function 
of the relative motion. Sec. [V| is devoted to discussions. 

II. FORMULATION 

A. Definition of the model 

The massive Schwinger model is two-dimensional QED with a massive fermion. It 
is not exactly solvable in contrast to the massless one The Lagrangian is given by 

C = -\f^ + $ [>f {id„ - eAJ -m]if>, (2.1) 

where = d^A u — d^A^. In two dimensions, the coupling constant e has mass dimension. 
It is therefore useful to measure all dimensionful quantities in units of e/ ypa. We hereafter 
set ej ^ = 1. Strong couplings correspond to small fermion masses. 

In the light-cone gauge (A + = 0), the only independent variable is ipn in the light-cone 
quantization. A~ and ipL are expressed in terms of iJjr as follows: 



'fe = ^. S , (2.3) 

with j + {x~) = a/2 : i/)' R (x~)il)ji(x~) :. We use the principal value prescription for (zcL) -1 
and (z<9_)~ 2 as in Refs. 



Eliminating A and ipL by using (|2.2| ) and (|2.3| ), one obtains the light-cone Hamiltonian 



p- = p- i p 

free ' i 

m 



free ' int > 

P 



/ cfe ) , (2.4) 



/ree 2^7- 
^ = 7T / ^-j+(x-)— — HaT) 



We expand ^ in terms of the creation and annihilation operators, 

1 poo fjh~^ r i , -| 

^«(«-) = 7TT77 / TTHfr fa^e - " *~ + ^(fc+K^I , (2.5) 



2V4 y (27r)\/fc+ 

where 6(fc + ) and d(k + ) satisfy the following anti-commutation relations, 

{b(k + ), = {d(k + ), = (27t)k + 8(k + - l + ) , (2.6) 

derived from |^(x~), ^(z/)} = {^/V^)8(x^ — y )- One may express P~ entirely in terms 
of b(k + ) and d(k + ) (and their Hermitian conjugates). We refer the reader to Ref. [|10| for 
the explicit form. 

We work in a truncated Fock space in which a state with total light-cone momentum 
P + — V is expressed as 

|V)p = |2)p + |4) 7 ,+ |6) 7J , (2.7) 
|2>„ = / . 1 2 8(h + k 2 - V)Mh, h)M\0) , 



^/(27r) 2 fcifc 2 

|4) P = - / LI -H=*(E - W^i, fc 2 ; fc 3 , fc 4 )6l444|0) , 
2 -A) i=1 ^(2tt)A; 4 i= i 

1 f^" 6 dk- 6 

|6) P = - / n ^^ 5 (E ^ - P)Mh, h, h; h, h, h)b\b\b\d\d\dl\Q) 
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where we use the abbreviated notations, b\ = b\ki), d\ = d){ki). We rescale momenta, ki — ► 
%i = ki/V and the wave functions, ip 2 (ki,k 2 ), i^iiki, k 2 ; k 3 , k 4 ), and ipQ^ki, k 2 , k 3 ; k 4 , k 5 , k 6 ) 
are replaced by ^2(^1,^2), 7 ?-1 V'40 i; i> x 2 ; x 3 , x 4 ), and P _2 ^ 6 (xi, x 2 , x 3 ; x 4 , x 5 , x 6 ), respec- 
tively. 

The wave functions ip 2 , ^4 and ip 6 must satisfy the following symmetry properties due 
to Fermi statistics, 

ip 4 (xi, x 2 ; x 3 , x 4 ) = -^4(^2, xi] x 3 , x 4 ) = -^4(^1, ^2; £4, £3) = ^4(^2, x 4 , x 3 ) , 
t/> 6 (zi, a?2, ^3; a?4, a?5, ^e) = -"06 (^2, x 3 ; x 4 , x 5 , x 6 ) = -ip 6 (xi, x 3 , x 2 ; x 4 , x 5 , x 6 ) (2.8) 

= -^6(^1,^2, £3! X5,X4,Z 6 ) = -06(^1,^2, ^3; X4, X 6 , X S ) etc. . 

If we require that this state has a definite property under charge conjugation transformation, 
we have further conditions on these wave functions, 

lfa(Xl,X2) = T'4 ) 2(x2,Xl) , 

ip 4 {xi, x 2 \ x 3 , x 4 ) = ±^4(23, x 4 ; x u x 2 ) , (2.9) 
ip 6 (x 1 ,x 2 , x 3 ; x 4 , x 5 , x 6 ) = T^e(x 4 , x 5 , x 6 ; x u x 2 , x 3 ) . 

The upper/lower sign in ( |2.9| ) corresponds to charge conjugation even/odd. 

The Einstein-Schrodinger equation M 2 \ip)-p = 2P~ P + \ip)-p leads to a set of complicated 
eigenvalue equations for the wave functions. It can be converted to a single matrix eigenvalue 
problem by expanding the wave functions in terms of basis functions, which we discuss in 
the next subsection. 

B. Basis functions 

It has been known that the wave function ip%(x, 1 — x) behaves as x 13 in the vicinity of 
x = ||, with P being the solution of the equation m 2 — 1 + nf3 cot(7r/3) = 0. By taking 
it into account, Mo and Perry concluded that a useful choice of the basis functions for the 
wave functions is given in terms of Jacobi polynomials, Pjf'^. In a previous paper fTJ], we 
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propose a simpler set of basis functions, essentially equivalent to that of Mo and Perry. We 
now propose another set of basis functions which leads to a drastic reduction of CPU time. 
We expand the wave functions as follows. 

N 2 

ip 2 (x, l-x) = ^2 a k F k (x, 1 - x) , 

k=0 

N 4 4 

i/ja(xi,x 2 ;x 3 , rr 4 ) = X&kG , k (xi,X2;x 3 ,x 4 ) , ^Xi = 1 , (2.10) 

k i=l 

N 6 6 

ip & {x u x 2 , x 3 ; x 4 , x 5 , x 6 ) = ckH k (x 1 ,x 2 , x 3 ; x 4 , x 5 , x 6 ) , ^ Xi = 1 , 

K i=l 

where we use the following basis functions: 

[x(l - x)f +k 

F k (x,l-x) = { (2.11) 
[x(l - x)f +k (2x - 1) , 

G k (x 1 ,x 2 ;x3,X4) = (xix 2 xzx A ) l3 (xix 2 ) kl (x' i x 4 ) k2 (xi - x 2 )(xs - x A )(xi + x 2 ) k:i , (2.12) 

H K (x 1 ,x 2 ,x 3 ;x 4 ,x 5 ,x 6 ) = {xxX-iX-iX^XfiY 

(x 1 x 2 ) kl (x l +x 2 ) k2 (x 1 -x 2 ) + (1,2,3 cyclic) 



x 



X 



(x 4 X5) fc3 (x 4 + a;5) fc4 (a;4-a;5) + (4,5,6 cyclic) (n + x 2 + x 3 ) k5 . (2.13) 



\kr, 



We abbreviate the upper limits of the sums. In reality, N 2 = 2 (Mi + 1), and it means that 
k in (gT[|) runs from to Mi. Similarly, iV 4 = (M 2 + 1) 2 (M 3 + 1), i.e., h, k 2 = 0, • • • , M 2 , 



fe 3 = 0, • • • , M 3 in ( pT2|) and iV 6 = (M 4 + 1) 2 (M 5 + 1) 2 (M 6 + 1), i.e., k 1: k 3 = 1, • • • , M 4 + 1, 
^4 = 0, • • • , M 5 , and k 5 = 0, • • • , M 6 in ( [2.13j ). The important point in choosing this set of 
basis functions is to reduce the number of the factors of the type (x± + x 2 ) k . This reduction 
allows us to express the basis functions in a simple way in the source code. For example, 
Gk(x\, x 2 \ x 3 , X4) may be written as 

G'k(x 1) x 2 ;x 3 ,x 4 ) ^B-l) 2 - 11 -' 2 ( 3 ] (x lX2 x 3 x 4 fx^x^x^x^ , (2.14) 

with iVi = k\ + i\ + j, N 2 = k x + 1 - i x + k 3 - j, N 3 = k 2 + i 2 , N A = k 2 + 1 - i 2 , and 
i = 2i\ + (We use a binary number for i.) Without using this new set of basis functions, 
six-body LFTD calculations would be much more heavy. 



We have explicitly separated the two-body basis functions into charge conjugation eigen- 
functions. But we have not done that for four-body and six-body basis functions because 
it makes the expressions so complicated that the drastic reduction of CPU time cannot be 
expected. We determine the charge conjugation property of an eigenstate by looking at the 
two-body state. From our experience, we know that it is a reliable way. 

With these expansions, the Einstein-Schrodinger equation becomes a (generalized) ma- 
trix eigenvalue problem, which can be solved numerically. Calculations of the matrix el- 
ements can be carried out analytically by using the formulas (and their generalizations) 
collected in an appendix of Ref . |J . 

III. NUMERICAL RESULTS 
A. Convergence 

First of all we have to see how many basis functions are enough to produce reliable 
results. We set m = 0.1 and gradually increase the number of basis functions. Fig. [1 
shows the lowest mass states in the calculation including only two-body basis functions. 
The lowest state is the meson state. It is charge conjugation odd. Its mass is 1.18160 at 
N2 = 10 (Mi = 4). The dashed line indicates the two- meson threshold. Note that the 
convergence is good enough for N2 = 10. 

As we increase the number of four-body basis functions, (keeping N 2 = 10) a state goes 
down below the two-meson threshold as shown in Fig. |2|. We regard this state as the two- 
meson bound state. On the other hand, the lowest state, the meson, is little affected by 
the inclusion of the four-body states. Its mass is 1.18103 at N4 = 80 (Mg = 3, M3 = 4), 
decreased only 0.05%. It is due to a negligibly small four-body component, 0.002%. Note 
that all the state above the two-meson threshold go down as N4 increases. We find that 
N4 = 80 is enough for the convergence for the lightest two states. 

We proceed to the six-body calculations, keeping N 2 = 10 and N± = 80 (N 2 + iV 4 = 90). 
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As seen in Fig. |3|, the convergence is quite good for the states below the three-meson threshold 
indicated by the dotted line. The three- meson threshold is given by the sum of the meson 
mass and the mass of the two-meson bound state. Again, the meson is not affected by 
the inclusion of six-body states. The two-meson bound state does not change either. The 
mass is 2.30980 at N 6 = 36 (M 4 = 2,M 5 = 1,M 6 = 0) which should be compared to 
2.31004 in the four-body calculations. The state just below the three-meson threshold can 
be regarded as the three-meson bound state, as we discuss shortly. It seems that Nq = 36 
(iV 2 + iV"4 + iV 6 = 126) is enough for the convergence for the lightest states. In the following, 
we restrict ourselves to this case. 

B. Two-meson bound state 

The second lightest state can be regarded as a bound state of two mesons. Its mass is 
2.30980 at m = 0.1. It has a 72.825% two-body component, 27.174% four-body component 
and 0.001% six-body component. The ratios change as the fermion mass changes. The 
smaller gets the fermion mass, the larger four-body component it has. For example, at m — 
0.01, it has 54.408% two-body component, 45.592% four-body component and 0.0004% six- 
body component. It is important to notice that it has a negligibly small six-body component. 
From this result we presume that it will not have any many-body components even if we 
could include higher Fock states. This state is charge conjugation even. 

C. Three-meson bound state 

We identify the three-meson bound state by the following criteria; (1) Its mass must be 
below the three-meson threshold. (2) It must have a large six-body component relative to 
the other states below the three-meson threshold, at least in the strong coupling region. (3) 
It must be charge conjugation odd. 

The first criterion is a trivial one, and is necessary for distinguishing it from three-meson 
scattering states. 
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The second criterion is based on the observation that the meson is almost completely in 
the valence state in the strong coupling region. As we discuss in detail in the next section, one 
may consider an (approximate) meson creation operator A\ Thus a three- meson state may 
be represented as ~ (v4t) 3 |0), which implies that the three-meson bound state has a large 
six-body component and negligibly small many-body components. Similarly a two-meson 
state has a negligibly small six-body component. 

The third criterion comes from the fact that a meson state is charge conjugation odd. 

We find such a state that satisfies all of these criteria. Its mass is 3.39181 at m = 0.1, well 
below the threshold, 3.49083. It has a 65.170% two-body component, 34.292% four-body 
components and 0.538% six-body component. All states near this have smaller six-body 
components, typically a few hundredths percent or less. For smaller fermion masses, the six- 
body component of the three-meson bound state become larger. For example, at m = 0.001, 
the two-, four-, and six-body components are 42.013%, 55.013%, and 2.974% respectively. 
It is charge conjugation odd. 

One might be surprised that it has a small six-body component, and might suspect that 
it is a two-meson state. In Sec. [V], we will argue that it cannot be regarded as a two-meson 
state. 

IV. WAVE FUNCTION FOR THE RELATIVE MOTION OF THE TWO-MESON 

BOUND STATE 

In the previous section, we utilize the information of the wave function to identify the 
three-meson bound state. It is an outstanding feature that we can get such information. 
In this section, we consider another example in which the information of the wave function 
is crucial. We introduce a meson creation operator to have a qualitative picture of the 
two-meson bound state by considering the wave function of the relative motion. 
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A. Meson operator 



Let us introduce an operator aJ(p), 

dk 



a\p) 



+ 



o (2ir)yjk{p-k) 
dk 



b\k)d\p~k) 



b\p + k)b(k) - d\p + k)d(k) 



>° (2tx) ^Jk{p + k) 
It is easy to show that it satisfies the following commutation relations 

a(p),a\q)\ = p5(p - q) , 

where a(p) is the Hermitian conjugate to a^(p) and annihilates the vacuum, 

a(p)|0) = 0. 

By using these operators, the Hamiltonian can be written in the following form, 



(4.1] 



(4.2) 



(4.3) 



An Jo k 2 



b\k)b{k) + d\k)d{k) 



1 r°° dp 

2 Jo 



+ - I ^a j (p)a(p) . 



(4.4) 



Note that in the massless case m = 0, the Hamiltonian is diagonal and the eigenstates are 
the Fock states of The operator a) is the creation operator of the meson in the (massless) 
Schwinger model, which is equivalent to a free massive scalar theory, 

r-P dk 



|meson(m = 0))-p = a\V)\0) 



-.b\k)d\V- fc)|0). 



(4.5) 



'o (2n)^/k(V-k) 

The meson is structureless in the sense that the wave function has no momentum dependence, 



ijj(k,V — k) = 1 for any k. Compare with ( |4.6| ) below. 

Once the fermion mass is introduced, the Fock states of are no longer the eigen- 
states of the Hamiltonian. The results in the previous section suggest, however, that one 
may introduce an approximate meson creation operator whose Fock states are approximate 
eigenstates of the Hamiltonian. We have seen that the wave function of the meson behaves 
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as if) ~ — x)} 13 and the higher Fock components are negligible in the strong cou- 
pling region. Taking into account these things, we propose the following approximate meson 
creation operator, 

fP dk 

A\p) = . ifj(k,p-k)bUk)dUp-k) (4.6) 

J o (2n)^k(p - k) 

POO (jh r 

+ / ip( P +k,k) b\p+k)b{k)-d ] {p+k)d{k) 

J o (27i)^k(p + k) 1 
where ip should be equivalent to ip 2 for the meson in the previous section and therefore had 
been known numerically, 

I meson (m ^ 0))?, ^ A\V)\Q) = / M,V - k)b\k)d\V - Jfc)|0). (4.7) 

(2n)Jk(V-k) 



o 

Compare with ( |2.7D . The shape of \tp\ 2 is shown in Fig. |5|. On the other hand, ip cannot 
be determined by looking at the meson state. But it affects two-meson states and can be 
determined by examining the two-meson bound state, at least in principle. 

It is hard to estimate the errors of the approximate operator A^(p), though the state 
( |4.7|) has been shown to be a fairly good approximation. For small fermion masses (small 
P), we expect that if) = 1 + 0(j3) and <p = 1 + 0(/3), and the errors are expected to be 0{(3) 
because A^{j>) reduces to a^(p) for if) = if = 1. 

B. Wave function of the relative motion 

Let us attempt to describe the two-meson bound state by introducing the wave function 
of the relative motion. Such a description is based on the assumption that it is a two-meson 
state, i.e., that the two mesons in the bound state would not be distorted too much. The 
assumption is justified a posteriori in the strong coupling region. 

Under this assumption, the two-meson bound state may be written as 

|two-meson) P = ^ f ^=5(Pi+P2 - VMp 1 ,p 2 )A\p 1 )A\p 2 )\0) (4.8) 

V2 JO y/piP2 

rV dk 

$ 2 (k,V-k)bUk)S(V-k)\0) 

o (2ir)yJk(V- k) 
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/ II / ' s(E h - V)Mki, H k 3 , h)b\bUld\\0) 
J o t=iJ(2ir)ki l= i 



2 

where we have substituted (|4.6|). The wave functions \l/ 2 and ^4 are expressed in terms of 
ijj and ip in the following way, 

* 2 (k h k 2 ) = -L I"' dq $(g, V - qMh, h - q)^(k 2 , h - q) (4.9) 

1 f k2 dq 
V2J0 J q ( V _ g ) 



*)(!/„./,,:/„./, /) - + fc 3, fe 2 + fc 4 ) ^ khk3) ^ hjh) ( 4.io) 

/ 2(A; 1 + fc 3 )(A; 2 + fc4) 



$(fcl+fc 4 ,fc 2 + fc 3 ) ,s 

+ , ^(fcl, fc 4 )^(fc 2 , fe 3 ) . 

/2CA:! + fc 4 )(A: 2 + fc 3 ) 



The wave function $ is that of the relative motion of the two mesons in the two-meson 
bound state. 

The meson operator ( f4.6| ) does not exactly satisfy the same commutation relations as 
( |4.2|) , but only approximately. Thus $(pi,p 2 ) does not need to be symmetric under the 
exchange of pi and p 2 . Nevertheless we regard it as being symmetric throughout this paper. 
Note that the ansatz (|4.8| ) is consistent with the charge conjugation symmetry, that is, 
^2(^1, k 2 ) is antisymmetric in k\ and k 2 , ^4(^1, k 2 ; k%, fc 4 ) is symmetric under the exchange 
of (ki,k 2 ) and (/c3,fc 4 ). 

It is interesting to note that this ansatz drastically reduces the degrees of freedom in 
the functional space. Due to the assumption that the meson wave function would not be 
distorted too much in the bound state, only one degree of freedom, i.e., the relative motion 
of the mesons, comes in. It seems natural to expand $(xi,x 2 ) (symmetric in x\ and x 2 ) as 

N 

Hx u x 2 ) = ]T Bi[ Xl x 2 ] l+l 2 , xi + x 2 = 1 , (4.11) 

1=0 

where Bi is a coefficient to be determined numerically. Taking into account the fact that 
ip(xi,x 2 ) is well approximated by ao(xix 2 )P, with being the normalization constant, 
a = [B{2(3 + 1,2/5 + l)] -1 / 2 (B is a Beta function), one may consider the following basis 
function expansions, 
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N 

V 3 (x,l-x) = Y,M x O- ~x)f +l (2x- 1) , (4.12) 
1=0 

N B a 2 

$?i{xi,X2, x 3 , x 4 ) = ^{xxxix^x^f (4-13) 
i=i V2 

x {-[(xi + x 3 )(x 2 + x 4 )f + [fa + x 4 )(x 2 + x 3 )] 1 } , 

where Ai is another coefficient to be determined numerically. Note that this set of basis 
functions is much simpler than the original one ( |2.12| ). 

We calculate the mass and the wave functions (i.e., the coefficients A\ and B{) of the two- 
meson bound state by using this set of basis functions with N = 7. The mass is calculated 
as 2.04180 at m = 0.01. This result is surprisingly good for this small set of basis functions. 
It is even better than the result of our full-set calculations, m = 2.05612. This is because of 
the factors like [(xi + Xs)(x2 + x 4 )] 1 in Q4.13Q , which are suitable for expressing the relative 



motion of the two mesons. It is therefore expected that this set of basis functions is good 
only for the two-meson bound state. 

The squared wave functions, |$(x, 1 — x)\ 2 , for various values of the fermion mass are 
shown in Fig. []. From this, we have an intuitive picture that the mesons are loosely bounded 
for small fermion masses, while they are close to each other for large fermion masses. This 
behavior has a simple physical interpretation: As is seen from Fig. |5], the meson is tightly 
bounded for small fermion masses. Therefore the very weak Van der Waals force between 
the two mesons causes a loosely bounded two-meson state. For large fermion masses, on the 
other hand, the meson has a broad shape. Therefore the charge distribution over a wide 
region keeps the two mesons close to each other by the Coulomb force. 

It is possible to check quantitatively how good the assumption (|4.8|) is. By inspection, 
we find that the wave function ip may be written as follows, 

<p(xux 2 ) = b (^y , (4.14) 

where the constant bo is very close to 1 for small fermion masses. This is consistent with 
the massless limit ( |4.1| ), in which /3 = 0, bo — 1, and if = 1. Given the form of if ( 4.14]) one 
may express \&2 i n terms of if), and if as 
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1 f xi dy 

^ 2 {x 1 ,x 2 ) = ^= / $(y,l-y)<p(xi,xi-y)ip(x 1 -y,x 2 ) (4.15) 

V2 J y (i - y ) 

--®(y, i - v)<p(x 2 , x 2 - y)ip(x u x 2 -y) . 



I / ' - dy 



V2h ^y(l-y) 

By substituting ip, ( |4.14|) , and Q4.11|) , we obtain the following expression for ^ 2 in terms of 



^2{XI,X 2 ) 



= V2 

x E (-i) 

fc=0 



N 



1 



B(l + m + 1,2/3+1) (4.16) 



1=0 m=0 

I + m — k 
k 



m 



(xix 2 ) 



where we used the formula, 



£*i*2-'= £(-i) 



n — m 



{xix 2 ) m , for xi + x 2 = 1 . 



(4.17) 



(=0 m=Q \ Til 

This should be compared with (|4.12 ). Fig. [7] shows the wave function ^> 2 calculated by the 
direct diagonalization and by eq. ( |4.16|) . (The coefficients B±, • ■ ■ , are obtained by the 
diagonalization. It is necessary to take into account the normalization condition, 



dx\<fr(x, 1 — x) 



(4.18) 



to obtain B . Although b still remains undetermined, we simply put fo = 1 for the compar- 
ison.) The agreement measures the validity of the concept of the "relative motion" of the 
two mesons inside the two-meson bound state. 



V. DISCUSSIONS 

By using a simpler set of basis functions, we have obtained the mass spectrum of the 
massive Schwinger model in the LFTD approximation. We have confirmed that the two- 
meson bound state has a negligibly small six-body component and found a candidate for 
the three-meson bound state. We emphasize that the information on the wave functions 
is very useful and is used for identifying the three-meson bound state. It is also used for 
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investigating the two-meson bound state. We introduce an (approximate) meson creation 
operator and the concept of the relative motion of the two mesons. This description gives 
an intuitive, qualitative picture of the two-meson bound state and motivates a very simple 
set of basis functions. 

The candidate for the three-meson bound state has a small six-body component compared 
to the two-body and four-body components, even in the strong coupling region. One might 
suspect that it is a four-body state, not a six-body state. Actually, the corresponding state 
appears below the threshold in the four-body calculations. Nevertheless, we think that it 
is the three-meson bound state for the following reasons. (1) In the massless theory, there 
exists only a free scalar particle in the physical spectrum. But because the creation operator 
of the meson, if expressed in terms of the creation and annihilation operators of the fermion 
and the antifermion as in ( |4.1|) , contains the annihilation operators, even a pure (free) three- 
meson state is not a pure six-body state. It is thus not so strange even if it has a small 
six-body component in the massive case. (2) The state is charge conjugation odd. A two- 
meson state should be charge conjugation even. It is hard to imagine that a two-meson 
state can be charge conjugation odd, if we rely on the description in terms of the meson. 
The description may be justified in the strong coupling region because the unperturbed 
(massless) theory has only the meson and the perturbation is small. It is natural to have a 
picture that the perturbation causes weak interactions between the mesons to form bound 
states. 

Unfortunately, we do not know why the six-body component of the three-meson bound 
state is so small. It is an outcome of complex non-perturbative effects. An analysis similar 



to that of Sec. |IVB| may reveal how the three-meson bound state looks like but will not 
explain the smallness of the six-body component. At this moment, we have to be content 
with showing that it is the three-meson bound state. 

In the strong coupling region, the two- and three-meson bound states appear above the 
threshold. We have the prejudice that despite our numerical results, they are, in reality, still 
bound. Probably they are just below the threshold in this region and approach the threshold 
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in the massless limit. We think that the reason why they do not appear to be bound is due 
to the limitation of our variational calculations and numerical errors. If one takes it for 
granted that they are really bound, one may estimate the errors in the calculations. 
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FIGURES 

FIG. 1. Two-body Tamm-Dancoff approximation. The lightest states are shown with the 
total number of basis functions. The fermion mass is m = 0.1. The lowest state is the meson. All 
the other states are "spurious." 

FIG. 2. Four-body Tamm-Dancoff approximation. The lightest states are shown with the 
total number of basis functions. The fermion mass is m = 0.1. The lowest state does not change 
at all. The second lowest state goes down below the two- meson threshold (the dashed line). It is 
the two-meson bound state. 

FIG. 3. Six-body Tamm-Dancoff approximation. The lightest states are shown with the total 
number of basis functions. The fermion mass is m = 0.1. The lowest two states do not change at 
all. The state shown by the line with triangle points is the candidate for the three-meson bound 
state. The three-meson threshold is indicated by the dotted line. 

FIG. 4. Fermion mass dependence of the mass eigenvalues. The dashed and dotted lines stand 
for the two-meson and three-meson thresholds respectively. 

FIG. 5. Squared wave functions for the meson, \ip(x, 1 — x)\ 2 , are shown for various values 
of the fermion mass. For small masses the wave function (in the momentum space) has little x 
dependence, implying that the meson is a compact object, while for large masses it is localized 
around x = 1/2, implying that the meson has a broad shape, and that the fermion and the 
antifermion are bound loosely. 

FIG. 6. Squared wave functions for the relative motion of the two mesons in the two-meson 
bound state, \&(x, 1 — x)| 2 , are shown for various values of the fermion mass. For small masses the 
wave function (in the momentum space) has a sharp peak at x = 1/2, implying that the mesons 
are bound loosely, while for large masses it has a round shape, implying that the mesons are very 
close to each other. 
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FIG. 7. Consistency check for the ansatz. The dashed line is the wave function *$>2(x, 1 — x) 
obtained by the direct diagonalization. The solid line is the wave function 1 — a;) constructed 

from <I> and \&4 by using the ansatz ( |4.8| ). The fermion mass is set m = 0.01. The agreement is 
quite good. It shows that our assumption for the two- meson bound state is quantitatively justified 
a posteriori. 
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